Theoretical models of ferromagnetic III-V semiconductors 
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Recent materials research has advanced the maximum ferromagnetic transition temperature in semi- 
conductors containing magnetic elements toward room temperature. Reaching this goal would make 
information technology applications of these materials likely. In this article we briefly review the 
status of work over the past five years which has attempted to achieve a theoretical understanding of 
these complex magnetic systems. The basic microscopic origins of ferromagnetism in the (III,Mn)V 
compounds that have the highest transition temperatures appear to be well understood, and efficient 
computation methods have been developed which are able to model their magnetic, transport, and 
optical properties. However many questions remain. 



It is hoped that the emerging field of semiconductor 
spintronics, which studies the controlled flow of charge 
and spin in a semiconductor, will lead to the develop- 
ment of new non-volatile, high-density, and high-speed 
information technologies. An important milestone in this 
field was the discovery five years ago of ferromagnetism in 
Mn-doped, p-type GaAs observed at temperatures in ex- 
cess of 100 Kelvin [1]. Ferromagnetism at room temper- 
ature with full participation of itinerant carriers would 
be a major breakthrough in semiconductor spintronics. 
With this aim, intensive material research is currently in 
progress on transition metal doped III-V semiconductors. 
In this brief review we present a snap shot of the theo- 
retical part of this endeavor, which has progressed hand 
in hand with experimental developments. 

A qualitative picture of the electronic structure of III-V 
diluted magnetic semiconductors (DMSs) was proposed 
by Dietl and coworkers [2,3]. Their model is based on 
the internal reference rule [6] which states that energy 
levels derived from the d-shell of the magnetic ion are 
constant across semiconductor compound families with 
properly aligned bands. The application of this rule to 
III-V materials is illustrated in Fig 1 for ionized magnetic 
impurities substituted on cation sites [2-5] . The position 
of the A2 level for Mn in GaAs, deep in the valence band, 
suggests that Mn in GaAs is 2+, that its d-shell is occu- 
pied by five electrons, and that incorporation of Mn in 
this and many other (III,Mn)V compounds will result a 
large density of valence band holes that can mediate fer- 
romagnetic coupling between the 5 = 5/2 Mn local mo- 
ments. For low carrier densities, these valence band holes 
will be bound to the Mn ion, leading to shallow accep- 
tor levels. This model of carrier-induced ferromagnetism 
is now fairly well established for (Ga,Mn)As as a result 
of electron paramagnetic resonance experiments [7], x- 
ray magnetic circular dichroism measurements [8,9], and 
magneto-transport data [10-13]. 



The acceptor impurity levels of Fe and Co [2,3,6,4] are 
unlikely to lead to high valence band hole concentrations 
in arsenides or antimonides, as shown in Fig. 1. The pos- 
sible coexistence of acceptor and neutral magnetic im- 
purities suggests that a ferromagnetic double-exchange 
mechanism can dominate in this case [14]. The scenario, 
however, has not yet been confirmed experimentally. For 
nitrides and phosphides the ionization energy of the ac- 
ceptor levels is even higher and the A2 Mn impurity level 
is not deep in the valence band, making it likely that 
d-charge fluctuations cannot be neglected. Experimen- 
tally, ferromagnetism with critical temperatures exceed- 
ing room-temperature has been claimed in (Ga,Mn)N 
[15,16] on the basis of hysteresis in magnetization mea- 
surements, however, the electronic configuration of Mn 
impurity in this compound has not yet been conclusively 
established, and these claims have not yet been widely 
confirmed. The observation of the p-type conduction 
achieved in some samples [17] seems to suggest that the 
simple picture of Ai Mn acceptor level according to Fig. 1 
may not apply in this case. 

Guided by the complex phenomenology very briefly 
discussed above and by available experimental data on 
III-V DMSs, several theoretical groups have performed 
detailed studies of these new ferromagnetic materials 
using both microscopic and effective Hamiltonian ap- 
proaches. The results of some of this theoretical research 
are reviewed in the following paragraphs. 

Ab-initio theories are an invaluable tool for study- 
ing the microscopic origins of ferromagnetism and for 
predicting electronic, magnetic, and structural ground- 
state properties of new compounds. The local density 
approximation (LDA) combined with disorder-averaging 
coherent-potential approximation (CPA) was used [18] to 
study GaAs and GaN DMSs doped with V, Cr, Mn, Fe, 
Co, and Ni. A ferromagnetic ground-state was predicted 
for V, Cr, and Mn doping while Fe, Co, and Ni impurities 
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were predicted to lead to a spin-glass state. The struc- 
tural properties of (Ga,Mn)As have been calculated in 
the dilute limit and for hypothetical zinc-blende MnAs 
using pseudoatomic-orbital-based LDA algorithms [19]. 
The calculated lattice constant of MnAs is similar to the 
lattice constant of GaAs. This result appears to suggest 
that the relatively large compressive strains observed in 
(Ga,Mn)As epilayers grown on GaAs are caused by in- 
trinsic defects rather than bv the substitutional Mn. 




FIG. 1. Approximate positions of transition-metal donor 
(D(0/+)) and acceptor (A(0/-)) levels in III-V compounds 
[3]. Acceptor state labeled by Ai and the deep A2 level were 
obtained from spin-resonance experiments in GaP:Mn [4] and 
from photoemission in (Ga,Mn)As [5], respectively. 

The electronic structure of (Ga,Mn)As supercells with 
Mn in the substitutional position, with As-antisites, 
or Mn-interstitial defects have been studied by vari- 
ous methods within the LDA [19,21,22]. Typical spin- 
polarizcd band and Mn 3d projected density-of-states re- 
sults are shown in Fig. 2 [21] for substitutional Mn. It 
has been conclusively established by the LDA calcula- 
tions that in the presence of substitutional Mn, the Fermi 
energy crosses the valence band near the T point, intro- 
ducing one hole per Mn. Clear theoretical evidence for 
As 4p-like band holes has been obtained in LDA+U cal- 
culations [20] and the Mn 3d density of states calculated 
by this method is in a good agreement with photoemis- 
sion data [23] . As-antisite and Mn-interstitial defects act 
as a double donor, according to the microscopic calcula- 
tions [19,21,22], and, therefore, are the likely sources of 
carrier compensation observed in experiment. In sum- 
mary, first-principles calculations support the picture of 
carrier induced ferromagnctism in some (III,Mn)V ferro- 
magnets, but detailed predictions depend on the approx- 
imation used for exchange and correlations and also per- 
haps on technical details of the calculation. This activity 
is on-going and further progress should be expected. 

Semi-phenomenological effective Hamiltonian theories 
in which the low-energy degrees of freedom are Mn ion 
S=5/2 spins and valence band holes, discussed in the fol- 
lowing paragraphs, have been developed to make quan- 



titative predictions for zero and finite temperature mag- 
netic, transport and optical properties of bulk samples 
[2,24] and heterostructures [25]. In (Ga,Mn)As, e.g., 
these models are partially justified by the negligible 
changes in the valence band structure of the host semi- 
conductor observed in angle-resolved photoemission ex- 
periments [26]. The key term in the effective Hamilto- 
nian description is a Kondo-likc exchange interaction, 
J p d J2i Sj • s(Rr), which arises microscopically from sp-d 
hybridization and strong d-orbital Coulomb interactions. 
Here S/ is the 5 = 5/2 local moment operator on Mn 
site / and s(r) is the multi-band envelope function hole 
spin-operator [27,28]. The simplest version of this model 
combines mean-field theory with a virtual crystal ap- 
proximation for the random distribution of Mn moments 
[2,29,24,30]. Illustrative results of the model are pre- 
sented in Fig. 3 for the critical temperature in (Ga,Mn)As 
and (In,Mn)As. These calculations and the model pre- 
dictions for the strain-engineered magnetic anisotropy, 
anomalous Hall conductivity [31], domain structure [32], 
and magnetic circular dichroism [27] are in good agree- 
ment with experiment. 
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FIG. 2. LDA Spin-polarized density of states of 
Gai4Mn2Asie crystal with substitutional Mn atoms [21]. 
Grey areas represent the contribution from Mn 3d states. En- 
ergy is measured from the Fermi level. 

The reliability of the standard mean-field approxima- 
tion has been examined by accounting for the Coulomb 
interactions among holes [24] , which enhance the critical 
temperature, and for correlations in Mn moment orien- 
tations which reduce the energetic cost of incompletely 
aligned spin configurations [33] and therefore lower the 
critical temperature. The calculations show that heavy- 
hole - light- hole mixing is responsible for a relatively large 
spin stiffness, resulting in only ~ 20 % reduction of the 
critical temperature due to spin-wave fluctuations in ar- 
senides and antimonides [34]. Moreover, the reduction 
is practically canceled by the enhancement of the criti- 
cal temperature due to hole-hole interactions (see inset 
of Fig. 3) which explains the success of the mean-field 
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theory in these ferromagnets [34]. 

Theoretical models beyond the virtual crystal approx- 
imation have been used to study the effects of disorder 
on transport properties, and magnetic properties in gen- 
eral, of the DMSs. The Boltzmann equation with Born 
approximation scattering rates have provided estimates 
of the anisotropic magnetoresistance effect (AMR) of or- 
der ~ 1 — 10 % [35], in good agreement with experi- 
ment [11,13]. AMR of a similar magnitude is predicted 
by Kubo formula for the ac-conductivity, as illustrated in 
Fig. 4. The key for understanding the transport and mag- 
netic anisotropy effects is a strong spin-orbit coupling in 
the host semiconductor valence-band. The absolute dc- 
conductivity estimates (see the left inset of Fig. 4) [35] 
are typically only 2-10 times larger than experimental 
values [11,12,36]. The discrepancy is likely partly due 
to inaccuracies in the model scattering amplitudes and 
multiple-scattering or localization effects omitted by this 
approximation. However experimental conductivities are 
still increasing as growth and annealing procedures are 
optimized so the discrepancy could be due mainly to un- 
intended extrinsic disorder. Kubo formula calculations 
of infrared conductivities, illustrated in the right inset of 
Fig. 4, of metallic (Ga,Mn)As DMSs demonstrate that 
transverse f-sum rule measurements can be used to ex- 
tract accurate values for the free carrier density [37]. The 
hole-fluid effective Hamiltonian theory discussed above 
uses no free parameters and is expected to be reliable in 
samples with high critical temperatures where the holes 
are metallic and their interaction with Mn and other im- 
purities is effectively screened. 
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FIG. 3. Main plot: Mean-field critical temperature calcu- 
lated using the effective Hamiltonian and virtual crystal ap- 
proximations for (Ga,Mn)As (triangles) and (In,Mn)As (cir- 
cles) DMSs as a function of the hole density [27,34]. Inset: 
Critical temperatures in (Ga,Mn)As obtained by including 
hole-hole exchange interaction (EX) and spin-wave fluctua- 
tions (EX+SW) is compared to the mean-field results (MF). 



-2 



4 


11 




l 1 l 


1 1 1 




- 1 ' l ' l 1 ' l 










1 3 


i 3a 0.94 Mn 0.06 A yJ 








- 




- B =-0.004 / - 








- 


"(32 












o 










- 


" — 1 
b 










- 


- 


r, i , i , i 




Ga 0.94 M(1 0.06 AS \ 




0,2 0,4 0,6 0,8 






hole density fnm" ) , 






=-0.004 \ 


- 








1 1,5 




\ - 








i | i i i | 










'A/) 1 

7 J 3 








/ • Ga 0.94 Mn 0.06 AS 


/ §.0,5 










♦ Ga 0.96 Mn 0.04 AS 

a Ga 98 Mn 02 As 




1 , 1 









400 800 
fira (meV) 

1 , 1 




-0,005 


0,005 


200 400 


600 800 





- 4 



2 

DC 



e nco(meV) 

FIG. 4. Left panel: dc anisotropic magnetoresistance 
(AMR) as a function of strain for 50% compensated (due 
to As-antisites) (Ga,Mn)As DMSs. AMR is defined as the 
relative difference between resistivities for in-plane magneti- 
zation oriented along the current direction and for magneti- 
zation perpendicular to the plane of the DMS epilayer. en 
is the relative difference between the substrate and the epi- 
layer lattice constants. Inset shows the absolute conductivity 
for magnetization along the current direction. Right panel: 
ac-AMR as a function of frequency. Inset shows the absolute 
ac-conductivity. 

More detailed studies of disorder effects in DMSs have 
combined Kondo description of the spin interactions with 
Monte Carlo techniques applied to both metallic [38] 
(high carrier density) and insulating [39] (low carrier 
density) regimes, and with the CPA [40] or spin-wave 
theory [41] for the metallic samples. The CPA calcu- 
lations, which start from a tight-binding approximation 
band model and use RKKY theory of the effective Mn- 
Mn interaction, predict a factor of 3 enhancement of the 
critical temperature due to disorder for typical Mn and 
hole densities in (Ga,Mn)As DMSs [40]. Other studies 
using dynamical mean field theory [42] suggest similar 
trends with some additional features when considering 
the effective J p d coupling as a variable. On the other 
hand, the random Mn distribution and spatial fluctu- 
ations of the magnetization suppress the ferromagnetic 
transition temperature in Monte Carlo simulations [38]. 
The non-collinear nature of the ground state due to dis- 
order has been pointed out in the spin-wave theory cal- 
culations [41]. It was also suggested [43] that highly frus- 
trated magnetic correlations may result from the strong 
spin-orbit interaction in the valence band and the spatial 
disorder. 

The metal-insulator transition regime was explored by 
exact diagonalization techniques [44,45]. Consistent with 
experiment, the studies predict a metal-insulator tran- 
sition for 1% Mn doping in GaAs DMSs [44,45]. It 
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has been emphasized [44] that defect correlations play 
an important role in the transition and the calculations 
[44,45] have provided evidence in favor of the applica- 
bility of hole-fluid and impurity-band magnetic-polaron 
models on metallic and insulating sides of the transition 
respectively. Monte Carlo studies of the insulating phase 
have found spatially inhomogeneous magnetizations bel- 
low the ferromagnetic critical temperature and an en- 
hancement of the critical temperature due to disorder 
[39]. 
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